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The NASA technical monitor for the documentation was Mr. T. L. Galloway. 

The Aerophysics Research Corporation project leader was Mr. D. S. Hague. The 
GASP program has been used by a number of companies and universities through 
NASA contracted studies and is under continuing development. Prospective users 
should consult NASA's Ames Research Center regarding the latest details of 
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1.1 INTRODUCTION 


Over the past several years, NASA's Ames Research Center has developed 
the General Aviation Synthesis Program, GASP. This computer program performs 
tasks generally associated with aircraft preliminary design and allows am 
analyst the capability of performing parametric studies in a rapid rnamner. 

GASP emphasizes small fixed-wing aircraft employing propulsion systems varying 
from a single piston engine with fixed pitch propeller through twin turboprop/ 
turbofan powered business or transport type aircraft. The program may be 
operated from a computer terminal in either the "batch" or "interactive 
graphics" mode. 

The program is comprised of modules representing the various technical, 
disciplines integrated into a computational flow which ensures that the inter- 
acting effects of design variables are continuously accounted for in the 
aircraft sizing procedure. The model is a useful tool for comparing configu- 
rations, assessing aircraft performance and economics, performing tradeoff and 
sensitivity studies, and assessing the impact of advanced technologies on 
aircraft performance and economics. By utilizing the computer model the 
impact of various aircraft requirements and design factors may be studied in 
a systematic manner with benefits measured in terms of overall aircraft perfor- 
mance and economics. 

The GASP program has as its purpose the numerical specification of many 
aircraft design characteristics. Input quantities are general indicators of 
aircraft type, size, and performance, and the synthesis is extended to the 
point at which all of the Important aircraft characteristics have been analyzed 
quantitatively. The synthesis model and procedure together develop the 
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aircraft oonf igura tions In a manner useful in parametric analysis and also 
provide a useful step toward more detailed analytical and experimental studies. 

The synthesis program consists of a control module and several technology 
subnodules which perform the various independent studies required in the design 
of general aviation or small transport type aircraft. Each of the six technology 
modules shown in Figure 1.1.1 is composed of one or more computer subroutines, 
and the input to each module may be either the output of another module, or it may be 
input directly to the module. The integrated approach ensures that results 
contain the effects of design interactions among the various modules. For 
example, a change in wing loading affects wing area, tail size, lift, drag, 
propulsion system size, cruise attitude, structural weight, range and other 
parameters. Any particular net effect may be large or small; nevertheless it 
is determined numerically regardless of its magnitude. 

1.1.1 Discussion 

This section provides a brief description of the engineering methods used 
in the synthesis program. The descriptions are in the order shown in Figure 
1 . 1 . 1 . 


X. 1.1.1 Geometry. In this module, the dimensions of the aircraft compo- 
nents are calculated. Typical input parameters are the number of passengers, 
aspect ratio, taper ratio, sweep angles and thicknesses of wing and tail 
surfaces. The cabin is assumed to be of circular cross section, and tail 
surfaces are sized using trend equations derived for existing aircraft. Out- 
put of this module provides areas, lengths, angles, etc., which may be needed 
by other mbdules. 
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FIGURE 1.1.1 - GASP PROGRAM STRUCTURE 













1. 1.1. 2 Aerodynamics. Lift coefficient is determined as the sum of a 
term proportional to angle of attack, and a term due to high lift devices 
such as slots, flaps, etc. Lift curve slope computation includes ground 
effect and the effects of aspect ratio, Mach number and sweepback. Drag 
coefficient is the sum of profile drag, increments due to high lift devices, 
landing gear and compressibility, and the induced drag due to lift, including 
ground effect. Configuration geometry, flight conditions and type of high 
lift devices are input, while drag polars are output for the cruise, takeoff, 
and landing flight condition. 

1. 1.1. 3 Propulsion. Currently, turbojet, turbofan, turboprop, and 
reciprocating or rotating combustion engines can be simulated. Both engine 
size and performance sure determined. Both cruise and take-off requirements 
of the aircraft may be specified. The results also provide engine thrust 
and fuel flow at any flight condition using performance data for the specific 
engine of interest. 

1. 1.1. 4 Weight and Balance. Gross weight and payload are input, together 
with details regarding aircraft geometry and weight trend coefficients. The 
program has options for sizing tip tanks and locating the wing siich that the 
aircraft is in balance for the center of gravity travel of the aircraft. An 
acceptable value of static margin is input for this purpose. 

1. 1.1. 5 Mission Performance. The taxi, take-off, climb, cruise and 
landing segments of a mission are analyzed, and total range is computed. 

Options are available for calculating engine out and accelerate/stop distance, 
best rate of climb, high speed climb and other operating characteristics. When 
a specific range is required, the aircraft size is determined which provides 
this range within a specified tolerance. 



1. 1.1. 6 Economics. Both flyaway and operating costs are determined in 


this moduli. Flyaway cost is found by summing estimates of labor costs, material 
costs, and purchased equipment costs including overhead, tooling, sales, and 
profit for manufacturer and dealer. Operating costs include fuel, oil, inspec- 
tion, maintenance, storage, insurance, depreciation, and taxes, and the variable 
and fixed costs are combined to determine total operating costs as a function 
of annual utilization rates. 

A typical computational flow through the GASP program is illustrated in 
Figure 1.1.2. 


1.1.2 Documentation 

The six major submodules of the GASP program, as listed in Figure 1.1.1 
are of quite different lengths and levels of complexity. In addition, many 
subroutines are called by more than one other subroutine, so that it may be 
unclear, for example, whether it is a "propulsion" or a "performance" subroutine. 
The choice is usually made arbitrarily, for the sake of convenience alone. 

The seven volumes of the report are organized as shown in Figure 1.1.3. 

The GASP program is composed of 65 computer subroutines, 48 of which are 
documented in detail. Utility subroutines are listed in Figure 1.1.3 for 
completeness j however, they were not documented in detail but are described 
in Section 1.3. 

Each of the subsequent volumes is organized by first defining the "major" 
and "minor" subroutines of that section. The discussion is then directed at 
explaining how the subroutines interact, and how the computer logic is related 
to the purpose of each subprogram. Each significant equation of the subroutine 
is defined! and discussed, and this discussion may include comment as to the 
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FIGURE 1.1.2 TYPICAL GASP COMPUTATIONAL SEQUENCE 
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Volume 1 - Control and utility Programs 
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assumptiona needed for its derivation, or other relevant detail. Part 2 of 
each volume is the user's manual in which the input and output parameters of 
each subroutine are tabulated and defined alphabetically, including the units 
in which each is measured. A sample problem is also represented in terms of 
its numerical input and output. Finally, Part 3 of each volume is a program- 
mer's manual showing detailed flow charts for all subroutines in the volume. 

1.1.3 Utility Subroutines 

The GASP system includes a number of subroutines which can be termed 
"utility" subroutines. These are relatively brief programs which may be 
called by several other subroutines, and which typically perform a numerical 
function such as tabular lookup. These utility programs cure listed alphabe- 
tically below, and are described very briefly in terms of their significant 
input and output quantities. 

The utility programs that will be used by the GASP system depends on the 
propulsion option being exercised. The utility programs may be catalogued by 
propulsion option as follows: 

1. Turbofan or propeller option uses the following: 


BIV 

MAXMHW 


INTS 

OUTPUT 


ITRLN 

TP ALT 


ITRMHW 

If the turbofan engine data is input in 

tabular form, then the 

following are used 

in addition to those 

in (1) above: 

BISC 

STORES 


DTABX 

TABX 


MAPS 

TTABX 
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3. If a propeller type of propulsion system is used, then the 

following are used in addition to those mentioned in (1) above: 

BILINE MAXBND 

BIQUAD UNINT 

1. 1.3.1 BILINE (T, I, XI, YI, Z, K) - Linear Interpolation, One 
Independent Variable . — Tabular interpolation generates a numerical value for 
Z, corresponding to input values of XI and YI. The tabular data T (I) specifies 
the table number; T(I + 1) » 0, 1 or 3 denotes the order of the interpolation; 
T(I + 2) is the number of X values; T(I + 3) is the number of Y values; and 
T(I 4) are the values of X in ascending order. Output K denotes the number 
of interpolations performed. 

1. 1.3. 2 BIQUAD (T, I, XI, YI, 2, K) - Quadratic Interpolation, One 
Independent Variable . — This subroutine performs an interpolation over a four 
point interval, to maintain slope continuity. Table number T(I), T(I + 1) 
is the number of X values; T(I + 2) is the number of Y values, and 

T(I + 3) sure the values of X in ascending order. Output K measures the 
number of interpolations. 

1. 1.3. 3 BISC (Y, X, N, IL, IH, J) . — This subroutine determines the "low" 

and "high" integers IL and IH specifying the output values Y<IL) and 
Y(IH) which bracket the input number x. The dimension of Y is N, and 
output J is 0, 1, or 2 according to whether Y(l) ^ X^Y(N), X-^Y(l) or 

X Y (N) respectively . 

1. 1.3.4 BIV (Z, X, Y, AX, AY, AZl, NX, NY, NERR) - Linear Interpolation , 

Two Independent Variables . — If input data X and Y fall in the tabular 
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range AX (NX) and AY (NY), respectively, then NERR « 1. The input data AZl 
is given at NX * NY points, and the output is Z unless X or Y fall 
outside the associated tabular range (X AX(1), etc.) in which case 
NERR = 2. 

1. 1.3. 5 DTABX (XTAB, YTAB, ZTAB, X, Z, L) This is a function which calls 

subroutines BISC and TABX, and which is itself called by TTABX. Independent 
variables X, Z define the dependent variable DTABX, according to principles 
of Lagrange interpolation. 

1. 1.3. 6 INTS (T, M, L, E, B, C, HMA, HMI, BET, DERIV) . — A finite difference 
integrator, performed in double precision, of a system of M simultaneous 
first-order differential equations, which are defined in external subroutine 
DERTV. The non-zero components of T(100) are related to the state variables 

in DERIV. The other parameters in the calling sequence are input, and sure 
associated with the numerical aspects of integration (error magnitudes, step 
sizes, etc.). 

1. 1.3. 7 ITRLN (AX, AY, X, Y, N) ♦ — This subroutine returns a value for 
Y corresponding to an input quantity X. The input p aram eters for the N 
pairs AX (IP and AY (I), and AX (I) must increase nonotonically. If X is 
less than AX(1) or greater than AX(N), the subroutine extrapolates for Y(X). 

1. 1.3. 8 ITRMHW (ERROR, ERRMl, DRIVER, F, FF, JC, JX) , Newton-Rapshon 
Method in GASP . — This subroutine determines a zero to a function defined 
externally. Inputs are ERROR, the current (non-zero) value of the dependent 
variable; DRIVER, the current value of the independent variable; and F, a 
multiplier near unity. Outputs are ERRMl and DRIVER, the augmented 
values of the dependent and independent variables, and JC, the counter. 

FF and JX are not used. 
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1. 1.3. 9 MAPS . — This program Is called by program MAIN, and It calls 
subroutine STORE 3 three times to develop tables for thrust, fuel flow and 
airflow in the cruise configuration. The independent variables sire altitude, 

Mach number and turbine inlet temperature ratio. 

1.1.3.10 MAXBND ( PARAM , PRMMl, DRIVER, DMIN, DMAX, F, FF, KC, KX) . — 
Determines the maximum values of the dependent variable PARAM, and the asso- 
ciated independent variable DRIVER, subject to DMIN — DRIVER ^ DMAX. F and 
FF are input multipliers near unity in magnitude, and KC and KX are output 
counters; KX is initially zero, and is set to 1 when the maximum is deter- 
mined. 

1.1.3.11 MAXMHW (PARAM, PRMMl, DRIVER, F, FF, KC, KX) . — This subroutine 
determines the maximum of an input function Y (X) “ PARAM (DRIVER) , which is 
defined externally. F and FF cure input multipliers near unity, and KC is 
an output interaction counter, while KX changes from 0 to 1 when the maximum 
is determined. The previous value of Y(X) is PRMMl, and DRIVER is both 
input and output value of X. MAXBND is similar to MAXMHW except limits are 
placed on DRIVER. 

1.1.3.12 OUTPUT . — This subroutine begins with thirteen common block 
statements, and it includes 34 FORMAT statements. The subroutine is called 
by MAIN for the purpose of printing over 100 input and output figures related 
to geometry, weights, aerodynamics of the aircraft. 

1.1.3.13 ST0RE3 (NMAPS , NPTS, NLINE, AMAP, Z, X, Y, IREAD, IPRINT, ITAPE) . — 
This is called by MAPS, and it stores the dependent variable Y(144, NMAPS) 

and the two independent variables X(12, NMAPS) and Z(12, NMAPS). Other 
input quantities are NMAPS the number of maps, NPTS, the number of points 


1-1 

10 







i 




» 


i 


on a line of constant Z; NLINE, the number of lines of constant Z, and 
AMAP, the identifying parameter of a map. The last three integers are also 
input, and at least one must be nonzero for the program to read or write data. 

1.1.3.14 TABX (XTAB, YTAB, 0, L) . — This function is called by TTABX, and 
it acts as an interpolation subroutine. In effect, TABX is the value of the 
independent variable XTAB (2) for which Y is zero, and this function calls 
subroutine BISC, which brackets the X-value 0 satisfying XTAB ( I ) ^ O^XTAB(J). 

1.1.3.15 TP ALT (ALTZ , ALT, PO, FKALT, TO GO, XKV) ■ — This subroutine 

relates static pressure, temperature and gravity, kinematic viscosity (PO, 

TO, GO, XKV) to the altitude. ALTZ is geometric altitude, ft. and ALT is 

potential altitude, ft, while PO is measured in lb per sq in. , TO in 

2 

deg R, and GO in ft per sec per sec. XLV is returned in ft per sec units. 

If PO is input, ALTZ and ALT cure output, and vice versa. FKALT determines 
whether geometric or geopotential altitude is used. 

1.1.3.16 TTABX (NMAPS , NPTS, NLINE, Z, X, Y, ZPR, XPR, WPR, ZVAL) , 
Interpolation, Three Independent Variables . — This is another function which is 
a four-dimensional interpolator, where NPTS are the number of points on a 
line, NLINE the number of lines on a map, and NMAPS the number of maps. 

For a choice of map value WPR, X-value XPR and z-value ZPR, the function 
takes the value TTABX. The dimensions are X(12, J) , Z(12, J) and Y(144, J) , 
where Y is the dependent variable and J is the map number. Typical inputs 
are values of temperature ratio, Mach number and altitude, and output might 

be thrust, fuel flow or airflow. 

1.1.3.17 UNINT (N, XA, YA, X, Y, C) . — This subroutine performs a four- 
point interpolation to generate a smooth curve with continuous slope between 
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adjacent intervals. The number of input pairs is N, and YA(I) is monotonic 
from I to ft. No such restriction applies to YA(I). If the input X is 
less than XA(1) , then let Y * YA(1); if X is greater than XA(N) then L ■ 2 
and Y * YA(N). Otherwise, L - 0 and Y is calculated by interpolation. 

1.1.4 External Subroutines 

The GASP program is composed of over 60 subroutines some of which call 
as many as 8 or 10 other subroutines. The alphabetic listing of these sub- 
routines is given in Figure 1.1.4 where the programs indicated parenthetically 
may be called by the indicated subroutine. The volume in which each subroutine 
can be found is also indicated in this tabulation. 

The contents of each volume of the documentation sure listed symbolically 
in Figure 1.1.5, where the parenthetic numbers correspond to the subroutines 
listed in Figure 1.1.4. 
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PROGRAMS AND THEIR SUBROUTINES 


PROGRAM 

MAIN (AEROUT, CTAER, DLAND, ENGSZ, ENGWGT, FLAPS, GACOST, MAPS, 

OUTPUT, PERFRM, PNOYS, RGBAL, SIZE, WGHT) 

SUBROUTINES - TURBOFAN AND PROPELLER OPTIONS 

1. ACCEL (DRAG, ENGINE, TP ALT) 

2. AERO 

3. AEROUT (CLIFT, DRAG) 

4. APPFLP (FLAPS, ITRMHW) 

5. ASPEED (CTAER, ENGINE, ITRMHW, TP ALT) 

6. BISC 

7. BIV 

8. CLIFT ..... 

9. CLIMB (CLIFT, DRAG, ENGINE, TP ALT) 

10. CTAER (AERO, CLIFT, DRAG, TP ALT) 

11. DERIV (CLIFT, DRAG) 

12. DLAND (AERO, CLIFT, DRAG, ENGINE, TP ALT) 

13. DLOAD 

14. DRAG (ITRLN) 

15. DTABX (BISC, TABX) 

16. ENGDTT (TTABX) 

17-23 ENGDTl-7 (ITRLN, BIV) 

24. ENGINE (ENGDTT, ENGDTl-7 , WACDG, ITRMHW) 

25. ENGSZ (APPFLP, DRAG, ENGINE, ENGWGT, PERFRM, TP ALT , TURN) . . . 

26. ENGWGT (ENGINE, HOPWSZ, RCWSZ) 


VOLUME 
. I 

. VI 
. Ill 
. Ill 
. IV 
. VI 
. I 
. I 
. Ill 
. VI 
. Ill 
. VI 
. VI 
. V 
. Ill 
. I 
. IV 
. IV 
. IV 
. IV 
. V 
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FIGURE 1.1.4 PROGRAMS AND THEIR SUBROUTINES 


SUBROUTINES - TURBOFAN AND PROPELLER OPTIONS (Continued) 


27. FLAPS (ITRLN, ITRMHW, TP ALT) Ill 

28. GACOST (ASPEED, ENGINE, TP ALT) VII 

29. INTS (DERTV) v I 

30. ITRLN I 

31. ITRMHW I 

32. MAPS I 

33. MAXMHW I 

34. NACDG IV 

35. OUTPUT (CLIFT, TP ALT) I 

36. PERFRM (ACCEL, CLIMB, DLAND, TAKOFF, TAXI, TURN, XRANGE) VI 

37. RGBAL (AEROUT, CTAER, ENGSZ, ENGWGT, FLAPS, OUTPUT, PERFRM, 

SIZE, WGHT) VI 

38. SIZE (TP ALT) II 

39. STORE3 I 

40. TABX I 

41. TAIL (BIV, CLIFT, ENGINE, ITRLN, TP ALT) V 

42. TAKOFF (CLIFT, DERIV, DRAG, ENGINE, INTS, TP ALT) VI 

43. TAXI (ENGINE, TP ALT) VI 

44. TPALT I 

45. TTABX I 

46. TURN (DRAG, ENGINE, TPALT) VI 

47. WGHT (DLOAD, ENGSZ, ENGWGT, TAIL) V 

48. XRANGE (ASPEED, CTAER, ENGINE, ITRMHW, TPALT) VI 


1-1 

14 



i 


FIGURE 1.1.4 PROGRAMS AND THEIR SUBROUTINES 
ADDITIONAL AND REPLACEMENT SUBROUTINES - USED BY PRQPRrrgp qptTONS 


PROGRAM 


BILINE 


BIQUAD 


ENGDAT (COST, GEARBX, PERFM, WAIT, ZNOISE) 


ENGINE (MAXBND, MAXMHW, PWRPLT, TP ALT, TURBEG) 


ENGSZ (APPFLP , DRAG, ENGINE, ENGWGT, ITRMHW, PERFRM, TP ALT) 


IV 


GEARBX 


HOPWSZ (ITRLN) 


MAXBND 


PERFM (BIQUAD, UNINT) 


PNOYS (ASPEED, ENGINE, GEARBX, TP ALT, ZNENG) . 


PWRPLT (ITRLN) 


TCWSZ (BIV, ITRLN) 


TURBEG (BIV, ITRLN, ITRMHW) 


UNINT. 


WAIT 


ZNENG (UNINT) 


ZNOISE (BILINE) 


» 



FIGURE 1.1.5 CONTENTS OF EACH VOLUME 


VOLUME 

CONTENTS 

I 

Introduction *(MAIN, 6, 7, 15, 29-33, 35, 39, 
40, 44, 45, 49, 50, 63) 

II 

Geometry * (38) 

III 

Aerodynamics *(2, 3, 4, 8, 10, 14, 27) 

IV 

Propulsion *(16-25, 34, 52-62, 65, 66) 

V 

Weight and Balance *(13, 26, 41, 47, 64) 

VI 

Performance *(1, 5, 9, 11, 12, 36, 37, 42, 43, 
46, 48) 

VII 

Economics *(28, 51) 


* Parenthetic numbers refer to subroutine numbers of 
Figure 1.1.4 
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Under 
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1.2 MAIN PROGRAM USER'S MANUAL 


Program MAIN acts as the control program in the computer synthesis 
of general aviation aircraft . By calling 14 principal subroutines , MAIN 
effectively controls all the 66 subroutines which make up the GASP package, 
and it is the input to MAIN which specifies the aircraft being designed. 

The GASP computer program is intended to apply to a broad spectrum of aircraft 
types, and each aircraft design is specified by over 200 aircraft input param- 
eters and about 60 propeller input parameters, as tabulated under namelists 
INGASP and INPROP in the following pages. 

Many different vehicle sizing and performance options are available in 
GASP. The user may select certain options and bypus others according to his 
needs by inputting appropriate values for several indicator variables . For 
example, economic and/or noise calculations will be performed or bypassed accor- 
ding to the values input for TBO and KNOYS, respectively. Likewise, mission 
performance calculations may be terminated at the end of any segment accor- 
ding to the value of IFLY. Thus, one of the important functions of Program 
MAIN is to control the sequence in which the various subroutines are called. 

It is obviously required that the input data be physically consistent, 
and for this reason the units of each input parameter should be carefully 
noted. Errors in input data will often be apparent in the numerical results. 
More troublesome, however, are those errors which have smaller, but still 
significant effects on the resulting design, since these errors may not be 
suspected. 

Many of the input parameters are given default values , and these are 
indicated parenthetically following the definition of the parameter . All other 
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parameters must be input before the program will run. Many variables are used 
only when certain program options are selected, and thus they need not be input 
when these options are not used. For example, 24 variables are used only when 
the tail is sized in TAIL (stability and control analysis, LCWING=2) . Likewise, 
some input variables are required only when noise and cost options are exercised. 
A typical input stream to GASP is presented in Figure 1.2.1. The data input 
begins with a Title Card assigned by the analyst which also contains the integer 
IENGSC. Data input basically follows the format: (1) Title Card; (2) Addi- 
tional data read by MAPS and STORE 3 , if IENGSC is negative on the title card; 

(3) NAMELIST INGASP; and (4) NAMELIST INPROP. 

The variables for these data blocks are presented in alphabetical order in 
Figures 1.2.2 to 1.2.4. Breakdowns of the Namelist INGASP and INPROP inputs 
arranged by categories are presented in Figures 1.2.5 and 1.2.6. Format of the 
Title card is 

TITLE CARD 

COL 2-72 used for title 

COL 75-76 Engine cycle indicator (IENGSZ) 

= 0, propeller aircraft (default value) 

= 1, General Electric CJ-610 
= 2, Garrett TFE 731 
= 3, UACL JT-15D 
= 4, AFCO/Lycoming ALF 502 
= 5, General Electric CF-34 
= 6, General Electric TF-34 
- 7, General Electric T700/F1-QCGAT 
= -1, engine data input in tabular form 

If engine data input in tabular form, engine data follows the title card and is 
set up as described in section headed Engine Table , Figure 1.2.7. 
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A typical input stream to the GASP program has the appearance 


uw turboprop specified cnciic sis 

• INGASP 

«• 13500 . ICS- 45 CMS. DF-J . NTYE-6. ENCRU- 4000. IMCRU- 10000.. 
K161ITE-2. 

SAB-3 . VS- 16 . AS- 1 . . WAS- 18 . PAX-19 . PS-40. 865. 

AR-7 71. TCP- 15. TCT-.IS. 0LNC4- 9. SLH- 400. TP- .334 
APHT-3 35. ARVT-I 544. TCVT- 09. TCWT- 09. SU1V-.35. SUM-. 40. 

WIG- 334. EVEV- 5. ELRV-0 . OELP-7.. CATD-0 . ALPH.0-2 . 

VMLFSL-357 . KNAC-I. SAW- 299. 

ELOON-2 216. EL00T-2 515. MCK-1 80. ELPC-4 73. 

WBARVX- 123. VBARKX-I . 165. C0CLTH- 235. B0EL TV- 1.639. 

0ELCO- 00I59.KWC0- 12. 

Aas-- 66 - 32.-.I0.0 .10.20. 30. .40. 575. .70. 90.1.13. 

AC DC OR- 1 8.1 175.1 05.1 035. 1 .009. 1.0.1 .009. 1 .025. 1 . 10. 1 .225, 1.85.2. IS. 
CFOC- 27. 8TE0B- 55. RaHAX-1 280. JTLTTP-4. 

Oa«TE ■ 56. OFLPTO- 0 . OFLPLD-36 . OCOOTE-.I2S. 

OELTT- 1500. KHUX-500 . 

JENGSZ-3. IPAPT-3. PC CPU* 0. . XLQOE-3 6. 

UUNAC-2 0. 

SKFW- 35. SKFS- 0256. SKPEl- 3575. SKLC- 04984. OWPQCH-30.. 

SKK- 4250. SKB-98 5.SKT- .243.SK2- . 356.UWPAX-I0I . . 

SKWW-155 I. 

WFEX-1749 . WTUL-686 . WPU-2300.. 

LCVINC-0. 

KAIL-0 . ICPUS-I. 

PWCPTX- 992. DVI-9.0. OVR-O.. XLTHAX-1 350. HU8-.360. 

OElTvR-l 0. 

RSMX-600 . XLFHX-I 10. T0ELAT-2. . TIOLE-300 . NTC-3.4. 

NCAOE-I . TB0-3000 . CMV-0.. CCRV-0 . SPPH- ISO . t 
or-o . FCSF- 70. OB-40 . CLI AS- 1000 . 

UCSENC-100 . ALP-50. i 

0VI-8 0. DVR. 7.0. OO.TVR-2.6. I 

ia«-3. vane -i 40 . 

KNAC-2. WENC-3S8 . WNAC-304 5 . 08ARN-2 91 .ELM-1 1 .04. 

KNAC-I. XLQOE-3 80. UWNAC -2 . 037 , SWSLS- 3143. SKPEl- .3575. FPTL-O.. 
OWACM- .35. CPALT-10000.. 

KM) 

6IKR0P 
HTTP- IS. 

WPP0PI-I5I.S. 

8L-3 .AT- 114.. 

at- s. 

XM1AX-41730 . 

KOOECP-4. 0PR0P-8.S. TSP0HX-I000.. KOOETM-6. 

FT.-1 0. 

KOOECR-7. MWSLS-640 . 

JSlZE-3. ANC(MP-.I2. 

WKPFAC-.863. 

KNOTS-l. MCTS-IOOO . 01 ST- 1000. 

too 


FIGURE 1.2.1 

TYPICAL INPUT STREAM 


ai-85. 

CP-04/92. 


OIST-IOOO. I DATE-1970, 


p °o*i?4 


Kg? 


This is card image of input deck: for propeller configurations both namelist 
"ingasp" and "inprop" are required; for turbofan configurations only namelist 
"ingasp" is required. 

Three examples illustrating use of the GASP program are presented in 
Appendices A, B f and C as follows: Appendix A - Turboprop Powered Design, 

Fixed Engine Size; Appendix B - Two-Place Trainer with Fixed Pitch Propeller, 
and Appendix C - Turbofan Design Using Scaled TFE-731 Engine. 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


VARIABLE 

ACDCDR 

ACLS 

ALPHLO 

ALTFLP 

ALTLND 

ALR 

AR 

ARHT 

ARVT 

ARVTE 

AS 

ATMXQC 

BENGOB 

BMLOD 

BOELTV 

BTEOB 

CATD 

CCRW 


DESCRIPTION 

normalized wing profile drag values in drag table 
(if KWCD ? 0 ) 

array of C L values in wing profile drag table (if 
KWCD + 0) 

zero lift angle of attack, deg 

altitude during takeoff and landing for Reynold's 
number calculation, ft (0.) 

altitude of landing field, ft (0.) 

manhour labor rate $ per hr (3.40) 

wing aspect ratio 

aspect ratio of horizontal tail 

aspect ratio of vertical tail 

effective aspect ratio of vertical tail (numerical 
function of ARVT and SAH) , if LCWING = 2 

number of aisles 

maximum tip tank length/wing tip chord (3.16 if KTIPX « 1) 

fraction of flap-free wing span due to engines (0.) 

length to diameter ratio of tail boon (14.5 if KC0NFG«1) 

wing span/vertical tail moment arm (if VBARVX input) 

flap span to wing span ratio (.75) 

0, normal design structural category, FAR Part 23 

1, utility design structural category, FAR Part 23 

2, aerobatic design structural category FAR Part 23 

3, transport design structural category FAR Part 25 

annual cost of crew, $(0.) 


CFOC 


flap chord to wing chord ratio (.3) 




FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


VARIABLE 


DESCRIPTION 


CHALF 


two-dimensional variation with angle of attack of ele- 
vator hinge moment coefficient (function of RH) if 
LOWING = 2 


CHDEL 


two-dimensional variation with elevator deflection of 
elevator hinge moment coefficient (function of RH) if 
LOWING » 2 


CINP 


cost of annual inspection, $ (1500.) 

fuselage form factor (numerical function of fuselage 
fineness ratio) 


CKHT 


horizontal tail form factor (numerical function of TCHT 
and SAH) 


nacelle form factor (numerical function of nacelle fine- 
ness ratio) 


CKTP 


CKVT 


CLEOC 


CLIAB 


CLTLMT 


CMFLPL 


tip tank form factor (numerical function of tip tan k 
fineness ratio) 

vertical tail form factor (numerical function of TCVT) 

wing form factor (numerical function of TCR and TCT) 

leading edge device chord/wing chord (0.) 

cost of liability insurance, § (215.) 

limiting C L in turn, if JTRSZ =» 1 (1.0) 

increment to fixed annual cost, $ (0.) 

wing C M about eg, landing flaps (function of DFLPLD) 
if LOWING » 2 


CMFLPT 


wing c about eg, takeoff flaps (function of DF1JPTO) 
if LOWING - 2 


pitching moment coefficient of all engines about eg at 
landing (0.) if LOWING ■ 2 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 







FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


VARIABLE 


DESCRIPTION 


DELTEO 


DELTT 


DELWFC 


DELTVR 


DELWST 


DEMAX 


DFLPTO 


DFLPLD 


DLMC4 


DLSWSH 


DRMAX 


DWPQCH 

DWPQCV 


EGMRGN 


ELINC 


optimum trailing edge flap deflection angle, deg. 

(default function of JFLAP) 

time spent taxiing before takeoff and after landing, hrs. 

incremental control group weight, lb. (0.) 

estimate of time required to rotate aircraft during 
takeoff, sec (3.5) 

incremental structural weight, lb. (0.) 

maximum up elevator deflection, deg (-25.), if LCWING=2 

takeoff flap deflection, deg 

landing flap deflection, deg 

sweep of wing quarter chord, deg 

increment in wetted area/wing area (0.) 

maximum rudder deflection, deg (25.0) if LOWING - 2 

increment of engine failure decision speed above stall, 
kts (5.) 

increment of takeoff rotation speed above engine failure 
decision speed, kts (5.) 

horizontal tail quarter chord sweep, deg, if LOWING/ 0 

vertical tail quarter chord sweep, deg, if lcwing/0 

aircraft depreciation period, year (8.) 

engine eg relative to leading edge of mac, for wing- 
mounted engines; fraction mac, positive aft (0.), if 
LCWING/0 

distance from leading edge of vertical tail to leading 
edge horizontal tail on line of intersection of vertical 
tail and horizontal tail, ft, if LOWING / 0 

nacelle length, KNAC - 2, ft 

length to diameter ratio of nose cone of fuselage (2.0) 




• FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 



VARIABLE 

DESCRIPTION 

• 

ELODT 

length to diameter ratio of tail cone of fuselage (3.2) 


ELPC 

length of pilot compartment, ft (4.44) 


ELRW 

length of pylon attachment, for fuselage mounted engines 

• 



EMCRU 

design cruise Mach number 


EMTURN 

turn Mach number, if JTRSZ = 1 


ENP 

number of engines 

ft. 

EYET 

horizontal tail incidence angle# deg (0.) if LOWING ** 2 


EYEW 

wing incidence to fuselage horizontal reference deg. 


FACWl 

change in gross weight to start range iteration 

• 


(default function of gross weight and range) 


FCSF 

fuel cost, $ per gal (.51) 

- 

FLAPN 

number of flap segments per wing panel (1.) 

ft 

FPYL 

factor for turbofan engine pylon weight (.7) if NTYE=7 
and KNAC/2 


FRESF* 

required reserve fuel; ^^10, fraction of 45 min;]>10, 
lb fuel (1.0) 

» 

GRFE 

landing gear flat plate area, sq'ftj (function of gross 
weight) 


HAPP 

landing obstacle height, ft (50.) 


HBTP 

turbofan engine face hub/tip ratio, if NTYE-7 and 

V 


KNAC /2 


HCK 

i 

mean fuselage cabin diameter minus mean fuselage nose 
diameter, ft (2.47) 

ft 

HIR 

hull insurance rate; insurance cost/ aircraft price (.02) 


HNCRU 

design cruise altitude, ft 

ft 

HOO 

altitude at start of mission, ft (0.) 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


I 


VARIABLE 


DESCRIPTION 


HPORT 


HTMAX 


HTURN 


HWING 


I CRUS 


IGEAR 


I PART 


I SWING 


takeoff altitude, when JENGSZ=1 or 2, ft (0.) 
hours between annual inspection (100.) 
wing height above ground during ground run, ft (3.) 
terminal altitude for takeoff segment, ft. (500.) 
altitude of turn, ft, if JTRSZ * 1 

0, low wing position on fuselage if LCWING = 2 

1, high wing position on fuselage if LCWING = 2 

1, climb at maximum rate of climb (default) 

2, climb at maximum allowable operating speed 

3, climb at input EAS 

0, cruise at EMCRU (default) for cost and range 
calculation 

1, cruise at normal power for cost and range 
calculation 

2, cruise for best specific range for cost and range 
calculation 

1, compute full mission (default) 

2, compute mission through takeoff segment only 

3, compute mission through climb segment only 

4, compute landing performance only 

type of landing gear: 

0, retractable (default) 

1, fixed gear 

1, FAR Part 25 Turbine (def atilt) propulsion Sizing 
requirements 

3, FAR, Part 23, General Aviation propulsion sizing 
requirements 

0, keep wing loading fixed during range balance 
(default) 

1, keep wing area fixed during range balance 

0, landing weight » gross weight (default) 

1, landing weight * weight at end of mission 

2, landing weight = fraction of gross weight 



FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 



VARIABLE 

DESCRIPTION 

ft 


0, size engine for cruise only 
l r size for cruise and takeoff 


JENGSZ 

2, size for cruise and takeoff and climb requirement 

3, size for cruise and climb requirement 

4, engine thrust specified; input KNAC = 2 , ELN, DBARN, 

ft 


WENG, WNAC , if NTYE = 7, only 

1, plain flap 

2, split flap 

3, single slotted flap (default) 

• 

JFLTYP 

4, double slotted flap 

JTRSZ 

5, triple slotted flap 

6, Fowler flap 

7, double slotted Fowler flap 

0, no turn (default ) (available only if NTYE-7) 


l r turn sizing option (available only if NTYE=7) 


KCONFG 

type of fuselage tail cone: 

0, conventional cone (default) 

1, tail boom support 

0, nacelle drag computed as penalty to. engine 



performance (turbo fans only) 


KNAC 

1, nacelle drag part of aircraft drag; nacelle sized 
by engine 


• 

2, same as 1, except nacelle size input DRARN, ELN 


EODETO * 

engine power setting during takeoff segment if NTYE=7 


KODECL * 

engine power setting during climb segment if NTYE- 7 


KODETR * 

engine power setting during turn segment if NTYE = 7 


KODEAC * 

engine power setting during acceleration segment if 
NTYE * 7 

* These variable are set to 5, 6, 7 where 


1 

5 = maximum power (default) 

6 = maximum continuous power 



7 = maximum climb power 


KPLOT 

0, no plotting (default) 

1 , aerodynamic data plotted 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 





FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 






FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


VARIABLE 

DESCRIPTION I 

SAH 

horizontal tail location on vertical tail: 1 

0, low tail I 

1, T-tail I 

SCFAC 

shift in divergence Mach number due to supercritical I 

design (0.) 1 

SINKTD 

landing touchdown sink rate, ft per sec (3.0) | 

SKB 

weight trend coefficient of fuselage (136.) 1 

SKCC 

weight trend coefficient of cockpit controls (11.) 1 

SKFS 

weight trend coefficient for fuel system (.0195) 1 

SKFT 

fraction of total theoretical tip tank volume used for 1 

fuel (.979) I 

SKFW 

weight trend coefficient of fixed wing controls (.404) 1 

SKLG 

weight trend coefficient of landing gear, fraction 1 

gross weight (.0318) 1 

SKMG 

weight trend coefficient main gear, fraction of landing I 
gear (.80) 1 

SKPEI. 

weight trend coefficient of engine installation, fraction 1 
dry engine (.135) 1 

SKPES 

weight trend coefficient of engine nacelle, fraction I 

dry engine (.338) if KNAC ? 2 1 

SKSAS 

weight of stability augmentation system, lb (0.) 1 

SKTL 

factor on tail weight for arresting hook (1.) 1 

SKWF 

fraction of total theoretical wing volume used for wing 1 
fuel (.430) 1 

SKWTP 

tip tank weight trend coefficient, lb per sq ft (1.89) | 

SKWW 

weight trend coefficient of wing without high lift 1 

devices (133.4) 1 

SKY 

weight trend coefficient horizontal tail (.18) 1 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN (INGASP) 


VARIABLE 


DESCRIPTION 


TDELLD 


temperature increment above standard during landing, 
deg. F, (0.) 


TDELTO 


takeoff temperature above standard during engine sizing, 
JENGSZ=1 or 2 (0.) 


TDELTX 


THEMAX 


takeoff temperature above standard during mission, 
deg. F (0.) 

maximum allowable fuselage floor angle, deg (15.) 

input thrust for one engine, lb., if JENGSZ=4 and 
NTYE=7 


TIDLE 


idle thrust for propeller configurations, lb., if 
NTYE^7 (0.) 


TROTID 


vertical position of thrust line relative to eg, positive 
for thrust below eg, ft (0.), if LOWING = 2 

property tax rate; tax/value (0.) 

ratio of reverse thrust to idle thrust during landing 

( 0 .) 


UCSENG 


UWNAC 


UWPAX 


VBARHX 


unit cost of engine, $ per lb thrust or $ per HP 
(default program calculates) 
coefficient of rolling friction (.02) 

nacelle weight/nacelle surface area; lb per sq ft, if 
KNAC*2 

weight per passenger, including baggage, lb (200.) 

horizontal tail volume coefficient (default function of 
fuselage length and diameter) 

vertical tail volume coefficient (default function of 
fuselage length and diameter) 

climb speed, EAS, kts (input only if ICLM-3) 

maximum structural design flight speed, mi per hr. 

ratio of allowable lift off speed to stall speed (1.1) 




FIGURE 1*2.2 INPUT - PROGRAM MAIN (INGASP) 



VARIABLE 

DESCRIPTION 

B 

VRATT 

ratio of landing approach speed to stall speed (1.3) 


VTDRAT 

ratio of touchdown speed to stall speed (1.15) 


WAS 

aisle width, in. 

B 

WCFLAP 

weight trend coefficient in flap weight equation 
(default function of JFLTYP) 


WENG 

dry weight of one engine, lb if KNAC = 2 


WFEX 

fixed equipment weight, lb (default function of PAX) 

Bo 

WFUL 

fixed useful load weight, lb 


WG 

initial gross weight, lb 


WGS 

wing loading, lb per sq ft 

B 

WLPCT 

ratio of landing weight to gross weight, if IWLD*2 


WNAC 

weight of one nacelle, lb if KNAC = 2 


WPLX 

design payload, lb (default function of PAX) 

B 

WPYLON 

weight of one pylon, lb, if KNAC =2 and TYPE 13 7 


ws 

seat width, in 

B 

WTMISN 

aircraft weight at start of mission, lb (default gross 
weight) 


WTRFAC 

weight during turn, % of gross, if JTRSZ * 1 (1.0) 


XLDGRQ 

required landing distance, ft (99999.) 

B 

XLFMAX 

maximum load factor during takeoff rotation (1.1) 


XLFMX 

i. 

landing flare load factor if < 4 or landing flare 
initiation height, ft if > 4. (1.2) 

B 

XLFTRN 

sustained turn load factor, if JTRSZ=*1 

XLQDE 

nacelle length to diameter ratio, KNAC»0 or 1 

B 

XTORQ 

required takeoff distance to clear 35 ft, input if 
JENGSZ=1 or 2 (99999.) 
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FIGURE 1.2.2 INPUT - PROGRAM MAIN { INGASP) 


VARIABLE 

YMG 


YP 


ZCG 


DESCRIPTION 


location of main gear on wing: 

0, on fuselage 

1 , at tip 

location of engines on wing: 

0, on fuselage 
1# at tip 

height above runway of eg at nose wheel lift off, ft, 
(function of HWING) if LOWING = 2 


j 







FIGURE 1.2.3 INPUT - PROGRAM MAIN (INPROP) 



VARIABLE 

DESCRIPTION 

• 

AF 

propeller or Q-FAN blade activity factor per blade 


ANCQHP 

nacelle area per horsepower (.12) 

• 

BL 

number of propeller or Q-E&N blades 


BLANG 

propeller blade angle at 3/4 rad., deg., only if speci- 
fying blade angle when NTYP = 1 


BMEP 

brake mean effective pressure, psi (0.) if NTYE > 10. 

ft 

CAMT 

initial production quantity of propellers to be used 
for costing (default function of propeller type) 


CLI 

propeller blade integrated design lift coefficient 

ft 

CTI 

initial estimate of propeller thrust coefficient (.2) 

ft 

DIST 

DPROP 

slant distance to observer for propeller noise , ft (1000.) 
propeller diameter, ft 


EMNOYS 

aircraft Mach number for noise calculation if KNOYS*0 

ft 

FT 

fraction of total propulsor thrust which is lost due to 
installation (0.) 


GR 

gear ratio, propeller rpm/engine rpm (1.) 


HCRIT 

critical altitude for turbocharger engines, ft (16000). 

ft 

HNOYS 

aircraft altitude for noise calculation, ft (1000.) if 
KNOYS 0 or 1 


HPM5LS 

maximum sea level static horsepower (0.) if KODECR-7 

ft 

HPQAB 

horsepower/bore area, piston engines, HP per sq in (2.6) 
if NTYE > 10. 

ft 

IDATE 

propeller weight technology level; 1970 or 1980# if 
NTYP > 10 


1-2 

17 



i 




FIGURE 1-2.3 INPUT - PROGRAM MAIN (INPROP) 


(CONTINUED) 


VARIABLE 


DESCRIPTION 


JSIZE 


KNOYS 


1, increase HP with constant propeller diameter 

2, increase both HP and diameter , keep disk loading 
constant (default) 

-1, no prop noise calculation (default) 

0, compute prop noise for aircraft at HNOYS and 
EMNOYS 

1, compute prop noise for aircraft at HNOYS at 
maximum level speed 


KODECR 


KODETH 


KSPCHG 


used during engine sizing 

for piston engines: 

m i # size engine and prop diameter at specified flight 
condition to maximize prop efficiency? engine 
operating point (PCRCR, PCPCR) specified 

for turboprop engines : 

= 1, engine being sized at a given flight condition? 

PCNCCR is input. T4 may be input T4STCR, otherwise 
T4/T2 = f (PCNCCR) 

for either piston or turboprop: 

= 2, size prop diameter at specified flight condition 
to maximize prop efficiency - engine size and 
operating point are fixed 

=3, 4, size engine at specified flight condition 

percent max engine power and prop size are fixed? 
KODECR = 3 , prop RPM not specified (iterate to 
max prop eff ) ; KODECR = 4, prop RPM is specified 
(no iteration involved) • 

s 7 , for horsepower and prop diameter input 

used during mission calculations 

3 5 or 6, find engine operating point (per cent max 
power) at specified flight condition for a fixed 
engine and prop size 

* 5, prop RPM not specified (iterate to minimize fuel 
flow) 

« 6, prop RPM specified (no iteration involved) 

0 f no turbocharger (naturally aspirated) (0) 

1 , turbocharged engine 





FIGURE 1.2.3 INPUT - PROGRAM MAIN (INPROP) 


(Continued) 


VARIABLE 

NCYL 

NTYP 


PCNCCL 

PCNCCR 

PCNCTO 

PCPCL 

PCPCR 

PCPTO 


DESCRIPTION 


number of cylinders, piston engines (4) if NTYE > 10 

1, fixed pitch propeller 

2, constant speed propeller 

3, constant speed full feathering propeller 

4, constant speed, full- feathering, deicing propeller 

5, constant speed, full feathering, deicing propeller 

with reverse 

6, Q-FAN propul sor 

11 to 16, same as 1 to 6, except Hamilton-Standard 
routines are used for estimating weight, cost, 
and noise 

per cent corrected rotor speed at climb for turboprop 

(1.0) if NTYE = 5 or 6 

per cent corrected rotor speed at cruise for turboprop 
(.96), if NTYE = 5 or 6 

per cent corrected rotor speed at takeoff for turboprop 

(1.0) , if NTYE = 5 or 6 

per cent maximum power in climb for reciprocating 
engines (1.0), if NTYE < 5 

per cent maximum power in cruise for reciprocating 
engines (.75) if NTYE < 5 

per cent maximum power at takeoff for reciprocating 
engines (1.), if NTYE < 5 




FIGURE 1.2.3 INPUT - PROGRAM MAIN (INPROP) 


(Continued) 


VARIABLE 

DESCRIPTION 

PCRCL 

per cent maximum rpm in climb for reciprocating engines 
d*) , if NTYE < 5 

PCRCR 

per cent maximum rpm in cruise for reciprocating engines 
(.907) , if NTYE < 5 

PCRTO 

per cent maximum rpm at takeoff for reciprocating engines 
(1.) , if NTYE < 5 

ROTN 

number of rotors , rotating combustion engines (2.) if 
NTYE = 14 

RWH 

ratio of width to height of piston engine (1.3) if 
NTYE > 10 

SKDIM 

dimension trend coefficient - engine cross-section (1.0) 

SKWGT 

weight trend coefficient - bare engine (1.0) 

T4STCL 

turboprop turbine inlet temperature at climb, deg R 
(Garrett TPE331 engine), if NTYE = 5 or 6 

T4STCR 

turboprop turbine inlet temperature at cruise, deg R 
(Garrett TYE331 engine) , if NTYE = 5 or 6. 

T4STTO 

turboprop turbine inlet temperature at takeoff, deg R 
(Garrett TPE331 engine) , if NTYE = 5 or 6 

TSPDMX 

maximum allowable propeller tip speed, ft per sec (900.) 

UCSPP 

unit cost of propulsor, $ per lb (default program calcu- 
lates) 

WKPFAC 

propeller weight adjustment factor (1.0) 

WPROPl 

weight of one propeller , lb . , if KNAC = 2 

XCLF 

propeller learning curve factor in costing for 1000 
units (1.02) 

XCLF1 

learning curve factor for single unit for propeller cost 
(3.2178) 

XCK70 

single unit propeller cost 1970 technology, $ per lb. , 
(default function of NTYP) 

XCKSO 

single unit propeller cost 1980 technology, $ per lb., 
(function of NTYP) 
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FIGURE 1.2.3 INPUT - PROGRAM MAIN (INPROP) 


(Continued) 


VARIABLE 


DESCRIPTION 


propeller counterweight factor (function of NTYP) 

coefficient in propeller or Q-FAN weight equation 
(function of NTYP) 

coefficient in Q-FAN shroud weight equation (function of 
NTYP) 

coefficient in gearbox weight equation (function of 
NTYP) 


XNMAX 


maximum engine speed, rpm 





FIGURE 1.2.4 OPTIONAL INPUT TO SUBROUTINES MAPS AND STORE3 


VARIABLE 

DESCRIPTION 

AMAP 

value of altitude 

IPRINT 

0, do not print input data 

1, print input data 

IREAD 

0, no data input 

1, read data from cards 

2, read data from Tape 11 

ITITL 

table title 

NLINE 

number of T4/T2 points 

NMAPS 

number of altitudes 

NPTS 

number of Mach number points 

SFNIDL 

idle specific thrust, lb per lb per sec 

T4MAX 

maximum turbine inlet temperature, deg R 

T4MC 

cruise turbine inlet temperature, deg R 

T4MCL 

maximum continuous or climb turbine inlet* temperature , 
deg R 

WAMAP 

SLS airflow of engine, lb per sec 

X (L, M) 

Mach number values 

Y(L, M) 

table values (thrust, fuel flow or corrected airflow) 
at altitude M 

Z(L, M) 

( 

T4/T2 values at altitude M 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF 
INGASP INPUT DATA 


GENERAL 
CONFIGURA- 
TION DATA 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 


WG 


Gross Weight (lb) 


WGS 

- 

Wing Loading (psf) 


PAX 


Number of Passengers (excluding pilot) 


ENP 

- 

Number of Engines 


IGEAR 

0 

= 0 - Retractable Gear; = 1 - Fixed Gear 


KCONFG 

0 

= 0 - Conventional Tail Cone 
= 1 - Boom Type Tail Support 


KTIPX 

0 

Tip Tank Indicator 

= 0 - No Tip Tanks; = 1 - Allows Tip Tanks 


NTYE* 

- 

Type of Engine Indicator 


KWRITE** 

- 

Print Control Parameter 


EMCRU 

- 

Design Cruise Mach Number 


HNCRU 

- 

Design Cruise Altitude 


*NTYE ■ 1 indicates reciprocating engine with carburetor. 

- 2 indicates reciprocating engine with fuel injection. 

- 3 indicates reciprocating engine geared with fuel injection. 

■ 4 indicates rotary combustion engine. 

* 5 indicates turboshaft engine. 

* 6 indicates turboprop engine. 

* 7 indicates turbojet or turbofan engine. 

=11, 12, 13 same as 1, 2, 3 except routine HOPWSZ used to compute 
engine geometry and weight. 

= 14 same as 4 except routine RCWSZ used to compute engine geometry and weight. 

■ 0 no propulsor. 

**KWRITE - 0 no print. 

■ 1 all write statements are printed. 

■ 2 selected summary write statements printed (normal output option). 

■ -1 selected summary write statements printed (abbreviated output option) . 

■ 9 additional write of propulsion performance (use for debugging). 


FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


GEOMETRY 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

FUSELAGE 

SAB 

- 

Seats abreast in fuselage 


ws 

- 

Seat width (inches) 


AS 

- 

Number of aisles 


WAS 

- 

Aisle width (inches) 


PS 

- 

Seat pitch (inches) 


ELPC 

4.44 

Length of pilot compartment (ft) 


HCK 

2.47 

Mean dia. cabin minus mean dia. nose (ft) 


ELODN 

2. 

Length/dia. ratio of fuselage nose section 


ELODT 

3.2 

Length/dia. ratio of tail cone 


BMLOD 

14.5 

Length/dia, ratio of boom (KCONFo *1) 

NACELLE 

KNAC* 

- 

Nacelle drag indicator 


ELN 

f(eng size) 

Nacelle length (KNA02) , ft 


DBARN 

f(eng size) 

Nacelle mean diameter (KNACs2) * ft 


ELRW 


Length of pylon attachment for fuselage 
mounted engines (ft) 

WING 

AR 

- 

Wing aspect ratio 


TCR 

- 

Wing root thickness/chord ratio 


TCT 

- 

Wing tip thickness/chord ratio 


SLM 

- 

Wing taper ratio 


DLMC4 

- 

Sweep of wing 1/4 chord (deg) 


EYEW 

- 

Wing incidence to horiz. reference (deg) 


*KNAC - 0 - nacelle drag accounted for in engine performance (only used with turbofans). 
" 1 - nacelle drag accounted for as an aerodynamic force; nacelle sized In 
engine routine. 

■ 2 - same as 1 except nacelle dimensions input in SIZE routine. 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA 


(Continued) 





HORIZ 

TAIL 

VBARHX 

f (geom) 



TCHT 

- 

• 


ARHT 

- 



SLMH 

- 



DWPQCH 

- 

• 


COELTH 

f (geom) 



SAH 



Horizontal tail volume coefficient 

Horizontal tail root thickness/chord ratio 

Aspect ratio of horizontal tail 

Taper ratio of horizontal tail 

Horizontal quarter chord sweep, deg 

Wing chord/horizontal tail arm 

Location of horizontal on vertical 
■ 0. - low tail; ■ 1 - T tail 


VERT 

TAIL 

VBARVX 

f (geom) 


TCVT 

- 


ARVT 

- 


SLMV 

- 


DWPQCV 

- 


BOELTV 

f (geom) 


Vertical tail volume coefficient 
Vertical tail root thickness/chord 
Aspect ratio of vertical tail 
Taper ratio of vertical tail 
Vertical tail quarter chord sweep, deg 
Wing span/vertical tail arm 


» 


I 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA 


AERO- 

DYNAMICS 


♦Form factor 
CKW - 1 

CKVT - 
CKHT - 

CKF - 1 

CKN - 1 
CKTP - : 


INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

CKW 

A 

Wing form factor 

CKF 

A 

Fuselage form factor 

CKN 

A 

Nacelle form factor 

CKVT 

A 

Vertical tail form factor 

CKHT 

A 

Horizontal tail form factor 

CKTP 

A 

Tip tank form factor 

ALPHLO 

- 

Angle of attack at C. « q 

Li 

DLSWSW 

0. 

Increment in wetted area/wlng area 

DELCD 

.0015 

Increment in C^ 

DELFE 

.25 

Increment In equiv. flat plate area of fuselage 

SCFAC 

0. 

0 - conventional drag divergence; 

> 0 - shift in Mp due to supercritical 

GRFE 

0. 

0 - correlated on gross weight; 

>0 - landing gear flat plate area (ft^) 

KWCD 

0 

Number of points in wing profile drag table 

ACLS 

- 

C^ values in wing profile drag table 

ACDCDR 

- 

Normalized wing profile drag values in 


wing profile drag table. 


defaults 

03 [2 + 4(t/c) w + 240(t/c)*] 

I + 4 ( 1 / 0 )^ + 240 ( 1 / 0 )^ 

1 + .lO(l-SAH) ] [2 + Mt/c)^ + 240(1/0)^] 

35 [1 + - -- ° + .0025(l/d)_] 

(1/Dj 

50 

+ ■ 60 v + .0025(l/d)_ 






FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


HIGH 

INPUT 



LIFT 

VARIABLE 

VALUE OF 


DEVICES 

NAME 

DEFAULT 

DESCRIPTION 


RCLMAX* 


Cl MAX ^asic wing at reference conditions 


ALTFLP 

0. 

Altitude for Reynolds number calc, ft 

FLAPS 

FLAPN 

i. 

Number of flap segments per wing panel 


WCFLAP 

f (JFLTYP) 

Coefficient in flap weight equation 


BENGOB 

0. 

Fraction of wing span without flaps due to 
wign mounted engines (0. - fuselage mounted) 


JFLTYP** 

3 

Flap type indicator 


DFLPTO 


Takeoff flap deflection, deg 


DFLPLD 


Landing flap deflection, deg 


CFOC 

.30 

Flap chord to wing chord ratio 


BTEOB 

.75 

Ratio of flap span/wing span 


DCLMTE 

f (JFLTYP) 

^ ' l MAX of ref * wing due to fla P s opt deflec 


DCDOTE 

f (JFLTYP) 

Ac d of ref. wing due to flaps at opt. deflec. 


DELTEO 

f (JFLTYP) 

Optimum flap deflection angle 

L.E. 

CLEOC 

0. 

L.E. device chord/wing chord ratio 

DEVICES 

DELLED 

0. 

Deflection of leading edge device 


DCLMLE 

.93 

ACl MAX of ref wing due to L * E ‘ devi ce at opt 


DELLEO 

45. 

Opt deflection angle for L.E. device (deg) 

> a __ TO. u J _ 1 a j.1. J , A < 


c/4 sweepback ■ Oo. Reynolds No. 
** Type of trailing edge devices: 
JFLTYP - 1, plain 
* 2, split 
■ 3 , single slotted 
« 4, double slotted 


6 x 10 6 


JFLTYP - 5, triple slotted 

■ 6, Fowler 

■ 7, double slotted Fowler 


This FLAPS routine is based on the methodology in the following reference: Sanders, 

Karl L. : "High Lift Devices, A Weight and Performance Tradeoff Methodology," 

Tech. Paper No. 761, The Society of Aeronautical Weight Engineers, Inc. May 
1969. 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF iNGASP INPUT DATA (Continued) 


INPUT 




PROPUL- VARIABLE 

VALUE OF 



SION NAME 

DEFAULT 

DESCRIPTION 

JENGSZ* 

- 

Engine sizing options 


IPART 

i 

1 - Part 25 turbine, 3 - Part 23 Gen Aviation 

PR 

i. 

Inlet pressure recovery factor 

THIN 

— 

Input thrust for one engine (lbs) (Input only 



if JENGSZ = 4) 


XTORQ 

99999. 

Required takeoff distance to 35 ft (input 



only if JENGSZ « 1 or 2) 


RWCRTX 

1.0 

Ratio of cruise wt/gross wt (used for eng. siz) 

RCCRU 


Required rate of climb @ cruise conditions 

HPORT 

0. 

Takeoff altitude (ft) | 

( 

Used only for 
engine sizing 

TDELTO 

0. 

Takeoff temp above std, 1 

(°F) 

when JENGSZ = 1 or 2 

SM1D 


Engine face Mach no. S.L. static ( Input only 




if KNAOO 

HBTP 

- 

Engine face hub/tip ratio 

or 1 

XLQDE 

- 

Nacelle length/diameter ratio | 

JTRSZ 

0 

0 = no turn, 1 = turn sizing option 

XLFTRN** 

- 

Turn load factor 


CLTLMT** 

1.0 

C T limit in turn 

i 



L 

1 Input only 

HTURN** 

- 

Altitude of turn 

) if JTRSZ-1 

EMTURN** 

- 

Turn Mach number 

( 

WTRFAC 

1.0 

weight during turn (% gross) or service ceiling 



weight 


ROCREQ 

50.0 

Engine out service ceiling rate of climb 

HSCREQ 

0. 

Engine out service ceiling required 


* JENGSZ ■ 0, size for cruise only 

■ 1, size for cruise and takeoff 

■ 2, size for cruise and takeoff 

and climb required 


■ 3, size for cruise and climb req. 

=4, engine thrust specified for turbo- 
fain aircraft (must use KNAC=2) . Must 
also input ELN, DBARN, WENG, WNAC 


** if turning performance is desired in mission profile, these variables must be 
input. Turning performance will be computed after climb segment. 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


INPUT 

PROPU- VARIABLE VALUE OF 

SION NAME DEFAULT DESCRIPTION 


Turbofan KODETO*** 5 Takeoff power indicator 

Version 

Only KODECL*** 5 Climb power indicator 

KQDETR*** 5 Turn power indicator 

KODEAC*** 5 Acceleration power indicator 


ft 




i 


ft 


*** if value “ 5, maximum power 

■ 6, maximum continuous power 

■ 7 , maximum climb power 


ft 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


WEIGHTS 


INPUT 

VARIABLE VALUE OF 

NAME DEFAULT DESCRIPTION 


SKPEI 

.135 

Wt coef engine instal (fraction of dry eng) 

SKLG 

.0318 

Wt coef landing gear (fraction of gross wt) 

SKMG 

.80 

Wt coef main gear (fraction of landing gear) 

SKPES* 

.338 

Wt coef eng nacelle (fraction of dry engine) 
0 for buried in fuselage 

SKY 

.180 

Wt coef horizontal tail 

SKZ 

.220 

Wt coef vertical tail 

SKTL 

1.0 

Factor on tail wt for arresting hook 

SKWW 

133.4 

Wt coef wing (excluding high lift devices) 

SKB 

136. 

Wt coef fuselage 

SKCC 

11. 

Wt coef cockpit controls 

SKFW . 

.404 

Wt coef fixed wing controls 

SKSAS 

0. 

Wt of stability augmentor system 

SKFS 

.0195 

Wt coef for fuel system 

SKWF 

.430 

Fraction of wing volume for wing fuel 

SKFT 

.979 

Fraction of theoretical tip tank volume 
for fuel 

SKWTP 

1.890 

o 

Tip tank wt coef (lb/surface area^ft ) 

LOWING 

0 

■ 0 - will not locate wing and balance 
aircraft; * 1 - balance aircraft 
• 2 - compute fwd and aft c.g. limits. 
Size tail based on stability and control. 

RELP 

0. 

Engine c.g. fraction of fuselage length (for 
fuselage mounted engines) 

EGMRGN** 

0. 

Engine c.g. in relation to L.E. of MAC 

(fraction of MAC) for wing mounted engines 

CPMRGN** 

.10 

Wing c.g. with respect to c/4 MAC 
(fraction of MAC) 
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FIGURE 1.2.5 . FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 



INPUT 

WEIGHTS VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 


STMRGN** 

0. 

Aircraft c.g. with respect to c/4 MAC 
(fraction of MAC) 


RELR 

.4 

c.g. of fuselage and contents (fraction of 
fuselage length) 


UWPAX 

200. 

Weight per passenger (UWPAX times PAX used 
maximum payload case) 


ATMXQC 

3.16 

Max tip tank length/wing tip chord 


LDCKMX 

8. 

Max 1/d of tip tank 


ELINC 

0. 

Distance between L.E. of V.T. and L.E. of 
H.T. on line of intersection of V.T. and 
H. T. (ft) 


WPLX 

f (PAX) 

Design payload (LB) 

• 

WFEX 

f (PAX) 

Fixed equipment weight (lbx) 


WFUL 

- 

Fixed useful load (includes crew) , lbs 


UWPAX 

200. 

Weight per passenger (UWPAX times PAX is 
used for maximum payload case) , lbs 


STRUT 

0 . 

Wing strut attachment point, fraction of 
semi-span (= 0, cantilever) 


VMLFSL 

- 

Maximum operating design flight speed (mph) 


CATD*** 

- 

Design category (structure) 

» 

DELP**** 

- 

Fuselage pressure differential (psi) 


YP 

YMG 


Location of engines on Wing. 0., on fuselage 
1 . , at tip 

Location of main gear on wing * 0 . - on fuselage 
and 1 . - at tip 


» 


* 

** 

*** 


Cooes frco ENGWGT routine if engine geometry computed there, otherwise default 
value 

Positive direction is aft; negative direction is forward. 

CATO ■ 0, normal (FAR 23) 


- 1, utility (FAR 23) 

- 2, aerobatic (FAR 23) 
■ 3, transport (FAR 25) 


Used to determine allowable load factors 
and design speeds 


**** if input DELP is not adequate to maintain an 8000 ft. cabin at cruise altitude 


the proper DELP will be computed in the program. If DELP is input as zero, 
it is assumed that the cabin is not pressurized. 


» 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


WEIGHTS 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

ENGINE 

WENG* 

WNAC* 

f(eng size) 
f(eng size) 

Dry weight of one engine, lb 
(Includes gearbox if geared) 

Wt of one nacelle, lb 


WPYLON* 

f(eng size) 

Wt of one pylon, lb 


SWSLS** 


Engine specific wt 

* lb/lb thrust for turbofan/jet 

* lb/HP for recip and turboprop 


UWNAC** 


Nacelle wt/nacelle surface area 
(lb/ft 2 ) 


FPYL** 

0.7 

Factor for pylon weight 


*Mu8t be Input if KNAC"2 (for non-zero weights); may be input for KNAC*0 
or 1 (no call to ENGWGT) . 

** Input only if KNAC-0 or KNAC-1. 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


PERFORMANCE 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 


HOO 

0. 

Altitude at start of mission, ft 


IFLY* 

1 

Partial mission indicator 


WTMISN 

WG 

Aircraft wt at start of mission, lb 


THEMAX 

15. 

Max allowable fuselage angle, deg 


UM 

.02 

Coefficient of rolling friction 


MUB 

.4 

Coefficient of braking friction 


HTG 

3.0 

Wing height above ground during ground run 

TAXI 

DELTT 

- 

Time spent to taxi before takeoff (hrs) 

TAKEOFF 

XLFMAX 

1.10 

Max load factor during takeoff rotation 


DELTVR 

3.5 

Guess on time required to rotate, sec 


DV1 

5.0 

Increment of decision speed above stall (kts) 


DVR 

5.0 

Increment of rotation speed above decision 


VRAT 

1.10 

Ratio of^ a^fiowable lift-off speed to stall r 


TDELTX 

0 

Increment in ambient temperature above 
standard day (°F) 


HTMAX 

500. 

Terminal altitude for takeoff segment, ft ACI 


NFAIL 

1 

- 0 - computes engine out and accel/stop dist 
• 1 - computes only all engine performance 

CLIMB 

ICLM 

1 

■ 1 - max rate of climb 

■ 2 - climb at max allowable speed 
- 3 - climb at input EAS (VCLMB) 


VCLMB 

0. 

Climb speed, EAS, kts 
(input only if ICLM«3) 


DELH 

1000. 

Altitude increment during climb 


*IFLY « 1 compute full mission 

" 2 compute mission through takeoff 
■ 3 compute mission through climb 
• 4 compute landing performance only 






FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


PERFOR- 

MANCE 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

CRUISE 

CRMACH 

EMCRU 

Cruise Mach number 


CRALT 

HNCRU 

Cruise altitude, ft 


ICRUS* 

0 

Cruise speed indicator 


FRESF 

1. 

Required reserve fuel 
< 10 = fraction of 45 min 
> 10 = lbs fuel 


RCRRQ 

0 . 

Required range or endurance 
= 0, no requirement 
<• 24, design endurance, hrs. 
>24, design range, n.mi. 


FACWI 

** 

Change in gross weight to start range iteratioi 


I SWING 

0 . 

= 0, hold wing loading fixed during range 
balance 

= 1, hold wing area fixed during range balance 

LANDING 

XLDGRQ 

99999. 

Required landing distance (ft) 


ALTLND 

0 . 

Altitude of landing field (ft) 


VRATT 

1.3 

Ratio of approach speed to stall speed 


RSMX 

1000. 

Maximum allowable rate of sink (fpm) 


TROTID 

0 . 

Ratio of reverse thrust to idle thrust 


HAPP 

50. 

Obstacle height (ft) 


SINKTD 

3. 

Touchdown sink rate (fps) 


XLFMX 

1.20 

Flare load factor (XLFMX<4) ; flare initiation 
height, ft (XLFMX £ 4) 


TDELAY 

1.0 

Delay for brake and reverse thrust application 
(seconds) 


IWLD 

0 

0, landing weight ■ gross weight 

1, landing weight » weight at end of mission 

2, landing weight = fraction of gross weight 


WLPCT 

- 

WLPCT, landing weight/gross weight ratio 
(IWLD = 2) 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 




» 


» 


INPUT 

PERFOR- VARIABLE VALUE OF 

MANCE NAME DEFAULT DESCRIPTION 


TDELLD 0. 


Temperature increment above std. (°F) 


VTDRAT 1.15 


Ratio of touchdown speed to stall speed 


» 


» 


* ICRUS ■ 0, cruise flown at input speed (EMCRU) 

>1, Cruise flown at speed at normal cruise power 

■ 2, Cruise flown at speed for best specific range 

** For gross weights below 5000 lbs or design ranges less them 800 n. mi., FACW1 ■ 
0.95. Otherwise FACWl * 0.75. 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


C0ST 


INPUT 


VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

NCADE 

0 

■ 0 no additional equipment cost 

■ 1 add equip cost a function of base cost 

CMV 

0 . 

Increment to hourly operating cost($) 

CCRW 

0 . 

Cost of crew ($) 

CMF 

0 . 

Increment to fixed annual co9t ($) 

HIR 

.02 

Hull Insurance rate (2/100) 

CL LAB 

215. 

Cost of liability insurance ($) 

PRV 

.20 

Aircraft residual value (2/100) 

DYR 

8. 

Years for depreciation (years) 

RI 

0 . 

Loan interest rate (2/100) 

TR 

0 . 

Property tax rate (2/100) 

CRW0H 

.50 

Crew overhead rate (2/100) 

CINP 

1500. 

Cost of annual inspection ($) 

HRI 

100. 

Hours between annual inspection (hrs) 

0HR 

5.5 

One engine overhaul cost ($/// T; $/HP) 

UCSENG 

f(NTYE) 

Unit cost of engine ($//*T; $/HP) 

UCSPP 

f (NTYP) 

Unit cost of propulsor ($///) (NTYP < 10) 

TB0 

- 

Time between overhaul (hrs) 

SRPM 

- 

Storage or tie-down rate ($/month) 

CP 

Routine 

Computes 

Aircraft price - if not input routine 
computes 

ALR 

3.40 

Manhour labor rate ($/hr) 

FCSF 

.51 

Fuel cost ($/gal) 








figure 1 . 2.5 functional listing of ingasp input data (continued) 
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STABILITY 
AND CONTROL 
TAIL SIZING 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 


LONGITUDINAL 

CMFLPL 

f <V 

Wing pitching moment coefficient about 
aircraft (landing flaps) 



CMFLPT 

£ <V 

Wing pitching moment coefficient about 
aircraft (takeoff flaps) 



CMPLD 

0 . 

Pitching moment coefficient about center 
of gravity due to all engines during landing 



STATIC 

.03 

Aircraft static margin, fraction of MAC 



CHALF 

f (RH) 

2-D variation of elevator hinge moment 
coefficient with angle of attack 



CHDEL 

f (RH) 

2-D variation of elevator hinge moment 
coefficient with elevator deflection. 



RH 

0.40 

Elevator chord/horizontal tail chord 



DEMAX 

-25. 

Maximum trailing-edge-up elevator 
deflection, deg «0 for T.E. up) 



EYET 

0 . 

Horizontal tail incidence angle relative 
to horizontal reference, deg 



ZCG 

ZAOf (HWING) 

Height of center of gravity above runway 
a nose wheel liftoff, ft 



TP 

0 . 

Vertical position of thrust line relative 
center of gravity, ft (> 0 for thrust 
below center of gravity) 



CXA 


Distance of main wheel contact point aft 
of MAC leading edge, fraction of MAC 

» 


DCMCLP 

0 for Jets 
f (T c ,DPROP) 

Propulsion stability term (d C /d C ) 
power, one engine m 



HWING 


Position of wing on fuselage 


■ 0, low wing 

■ 1, high wing 
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FIGURE 1.2.5 FUNCTIONAL LISTING OF INGASP INPUT DATA (Continued) 


STABILITY INPUT 


AND CONTROL 
TAIL SIZING 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

LONGITUDINAL 

TAUH 

f (RH) 

Elevator effectiveness 

DIRECTIONAL 

CNPAC 

f (WG, B) 

Required directional stability of aircraft 
C n^, per deg 


ARVTE 

f (ARVT, SAH) 

Vertical tail effective aspect ratio 


RV 

0.40 

Rudder chord/vertical tail chord 


RVMCS 

1.0 

Minimum control speed 

Stall speed (takeoff configuration) 


DRMAX 

25. 

Maximum rudder deflection, deg 


Rudder effectiveness 

da Vl/ d ® Rudder 


TAUV 


f (RV) 


FIGURE 1.2.6 FUNCTIONAL LISTING OF INPROP INPUT DATA 




INPUT 

VARIABLE 

VALUE OF 


m 

ENGINES 

NAME 

DEFAULT 

DESCRIPTION 

w 


KODECR* 

- 

Recip/turboprop engine cruise sizing option 



KODETH* 

- 

Recip/ turboprop eng throttling options 

• . 


XNMAX 

- 

Max engine speed, rpm 



GR 

i. 

Gear ratio = propeller spd/eng spd 



HPMSLS 

0. 

Max SLS horsepower; input if K0PECR-7 

• 


ANCQHP 

.12 

Nacelle area/horsepower (for NTYE 10) 



JSIZE 

2 

Engine sizing indicator, takeoff and climb: 
= 1 increase HP with no inc in prop diam; 

= 2 increase both power and prep diam but 





hold disk loading const (lll'MSI S/AIHSK) 

• 






* KODECR - used during engine sizing 

For piston engines: 

f KODECR = 1, size engine and prop diameter at specified flight condition 

to maximize prop efficiency; engine operating point 
(PCRCR, PCPCR) specified. 

For turboprop engines: 

KODECR = 1, engine being sized at a given flight condition; PCNCCR 

% is input. T4 may be input T4STCR, otherwise T4/T2 = f (PCNCCR). 

For either piston or turboprop: 

KODECR = 2, size prop diameter at specified flight condition to 

maximize prop efficiency - engine size and operating point 
are fixed. 

£ =3, 4 size engine at specified flight condition - per cent max 

engine power and prop size are fixed; KODECR = 3, prop 
RPM not specified (iterate to max prop eff ) ; KODECR = 4, 
prop RPM is specified (no iteration involved) . 

= 7, for horsepower and prop diameter input. 

• * KODETH - used during mission calculations 

KODETH = 5, 6 find engine operating point (% max power) at specified 
flight condition for a fixed engine and prop size; KoDETIi 
= 5 prop RPM not specified (iterate to minimize fuel flow; 
KODETH = 6 , prop RPM specified (no iteration involved) . 
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FIGURE 1.2.6 FUNCTIONAL LISTING OF INPROP INPUT DATA 


ENGINES 

INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 

RECIP 

PCPTO 

1. 

% power @ takeoff for recip engine 
(= POWER/HPMSL) 


PCRTO 

1 . 

% RPM @ takeoff for recip engine 
(= RPM takeoff /XNMAX) 


PCPCL 

1 . 

% power @ climb for recip engine 

(= POWER /HPMSLS) 

CL 


PCRCL 

1. 

% RPM @ climb for recip engine 

(= RPM /XNMAX) 

CL 


PCPCR 

.75 

% power @ cruise for recip engines 

(= POWER /HPMSLS) 

CR 


PCRCR 

.907 

% RPM @ cruise for recip engine (=RPM / 
XNMAX) CR 


KSPCHG 

0 

Supercharger indicator: 

= 0, naturally aspirated engine 
= 1, supercharged engine 


BMEP 

0. 

Brake mean effective pressure, psi 


HCRIT 

16000. 

Cricital altitude, ft (KSPCHG=1) 

TURBOPROP 

PCNCCR 

0.961 

% corrected rotor speed at cruise 
( turboshaft/prop) 


PCNCCL 

1.0 

% corrected rotor speed at climb (turbo- 
shaft/prop) 


PCNCTO 

1.0 

% corrected rotor speed at takeoff 
{ turbos ha f t/prop ) 


T4STCR 

0. 

Turbine inlet temperature at cruise, °R 
(turboshaft or turboprop) 


T4STCL** 

0. 

Turbine inlet temperature at climb, °R 
(turboshaft or turboprop) 


T4STTO** 

0. 

Turbine inlet temperature at takeoff, °R 
(turboshaft or turboprop) 


** If the default values (zero) are used, the program uses the limits specified 


in routine TURBEG for the Garrett TPE 331 turboprop. 
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FIGURE 1.2.6 FUNCTIONAL LISTING OF INPROP INPUT DATA 


9 

INPUT 

VARIABLE 

PROPELLER NAME 

VALUE OF 
DEFAULT 

DESCRIPTION 


NTYP* 

- 

Type of propulsor indicator 

• 

AF 

- 

Propeller blade activity factor/bladc 


DPROP 

- 

Propulsor diameter, ft 


BL 

- 

Number of propeller blades 

• 

CLI** 


Prop blade integrated design lift coefficient 


BLANG 

* 

Propeller blade angle @ r/R=.75 (deg) 
(this is only input if blade angle is 
specified for fixed pitch) 

9 

IDATE 

- 

Propeller tech level, 1970 or 1980 


TSPDMX 

900. 

Max propeller tip speed, ft/sec 

9 

FT 

0 . 

Thrust loss factor (fraction of total 
thrust: T = (1 - FT) T F(r _Q * FT =-1-0, 

Program computes FT. 


CTI 

.2 

Initial guess on thrust coefficient 
(propeller) 

9 

PCLER 

0.058 

Propeller tip - fuselage clearance, 
fraction of propeller diameter 


* NTYP = 1, fixed pitch propeller 

= 2, constant speed propeller 

= 3, constant speed, full feathering propeller 
= 4, constant speed, full feathering, de-ice propeller 
= 5, constant speed, full feathering, de-ice propeller, with reverse 
= 6, QFAN propulsor 

= 11, 12, 13, 14, 15, 16 - same as 1, 2, 3, 4, 5, 6, except Hamilton 
Standard routines used for propulsor weight, cost and n oisc. 

** Recommended value: CLI = .5 


t 
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FIGURE 1.2.6 FUNCTIONAL LISTING OF INPROP INPUT DATA 



INPUT 

VARIABLE 

NAME 

VALUE OF 
DEFAULT 

WEIGHTS 

XKl 

f (NTYP) 

(input only 
if KNAC=1) 

XK2 

f (NTYP) 


XK3 

f (NTYP) 


XCW 

f (NTYP) 


BNUM 

- 


AFT0T 

- 


SKWGT 

1.0 


SKDIM* 

1.0 


RWH 

1.3 


NCYL 

4 


R0TN 

2. 


HPQAB 

2.6 


WKPFAC 

1.0 


WPR0P1 

- 

COST 

SCLF1** 

3.2178 


XCIJ'** 

1.02 


XCK70** 

Computed 


SCK80** 

Computed 


CANT** 

Computed 

NOISE 

KN0YS*** 

-1 


I' 1ST 

1 000 . 


HN0YS 

1000. 


EMN0YS 

_ 


DESCRIPTION 

Coefficient in propulsor wt rqu;it ion 

Coefficient in propulsor shroud wt 

Coefficient in gearbox wt 

Accounts for propeller counterweights 

Number of blades for propulsor (QFAN) 

Total activity factor of QFAN 

Wt coefficient - bare engine wt 

Dimension coef - eng crosssoc t iona l dimension 

Ratio of width/height of piston engines 

Number of cylinders - piston engines 

Number of rotors - R/C engines 

HP/bore area - piston eng (HP/in") 

Propeller wt adjustment factor 

Wt of one propeller, lb (KNAC--2 only) 
(includes gearbox, if geared) 

Learning curve factor for single mi: it 
Learning curve factor for 1000 uni I s 


Single 

unit 

O.E.M. 

prop cost L-170, 

$/u 

Si nglo 

uni t 

0 . E . M . 

prop cost 

$/ii 


Initial quantity to be used 
Propeller noise indicator 
Slant distance to observer, M 
Aircraft altitude for noise rale., ft 
Aircraft Mach no. for noise rale. (KN0YN 0) 


^Diameter for rotary combustion engines; width for piston engines 
‘^Default values for these parameters are taken from NASA CR-2066, "Computer 

Program User’s Manual for Advanced General Aviation Propeller Study,” Muv l c J7 r > 
*** KN0YS = -I No noise calculation 

= 0 Compute noise for aircraft flying at UN0YS and KMN0YS 

= 1 Conput noise for aircraft flying at max level speed at IIN0YS 
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FIGURE 1.2.7 





INPUT FORMAT FOR ENGINE TABLE (IF 
TURBOFAN VERSION ONLY 

I ENGS 2 = -1) 


Card 1 





• 


IREAD 

Col. 

1-5 

right j u s t i f S e 



IPRINT 

Col. 

6-10 

right justific 



WAMAP 

Col. 

21-30 

left justified 

• 


T4MAX 

Col. 

31-40 

left justified 



T4MCL 

Col. 

41-50 

left justified 



T4MC 

Col. 

51-60 

left justified 

• 


SFNIDL 

Col. 

61-70 

left justified 


ft 






ft 


ft 


IREAD = 0 No data input 

- 1 Read data from cards 
= 2 Read data from Tape 11 

IPRINT =0 Do not print input data 
= 1 Print input data 

WAMAP = SLS Airflow of engine (]b/sec) 

T4MAX = Maximum turbine inlet temperature (°P*) 

T4MCL = Maximum cont. or climb TIT (°R) 

T4MC = Cruise TIT (°R) 

SFNIDL = Idle specific thrust (Ib/lb/sec) 

Card 2 

ITITL - Table Title 


Card 3 

NMAPS - Number of altitudes 


i 


i 
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Card 4 


Blank Card 


1 


Card 5 


NPTS - Number of Mach number points 
NLINE - Number of T4/T2 points 
AMAP - Value of Altitude 


Card 6 

X(l, 1) X(NPTS , 1) - Mach no. values 


Croup 

for 

each 

Altitude 


Card 7 

3(1, 1), Y (L, 1) Y (NPTS j 1) 


2(NLINE, 1), Y(NLINE, 1) 


2 * T4/T2 values Y 
(thrust, fuel 


= Table values I 

flow, or airflow) j 
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APPENDIX A 


TURBOPROP POWERED DESIGN, FIXED ENGINE SIZE 
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END OF ACCELERATION segment 
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DC OF ACCELERATION segment 
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1.3 PROGRAMMERS MANUAL FOR 


MAIN PROGRAM AND UTILITY SUBROUTINES 

This volume presents a description of the GASP Executive Program 
and the Utility Subroutines of GASP . 

1.3.1 MAIN Program 

The principal purposes of the MAIN program are the reading of input data 
required for the aircraft design, and the calling of the subroutines which 
carry out this design. The input data is read as a title card , NAMELIST/ 
INGASP/ and NAMELIST/INPROP/ , and tabular input if that propulsion option 
is selected which total about 220 and 50 parameters respectively. Many of 
these are given default values in the event that no numerical value is 
assigned in the NAMELIST format. These parameters vary from the fundamental 
(gross weight, cruise Mach number, etc.) to the detailed (takeoff rotation 
rate, seat width, etc.), and are listed alphabetically in Section 1.5. The 
subroutine structure of MAIN down to the first level arrayed by technology 
is presented in Figure 1.3.1. 

The main program calls one minor data reading subroutine (MAPS) and 
thirteen major subroutines which are normally called in the following order. 
Each subroutine may call other subroutines as indicated parenthetically: 

SIZE 

FLAPS 

DLAND (AERO, CLIFT, DRAG, ENGINE) 

CTAER (AERO, CLIFT, DRAG) 

ENGSZ (APPFLP , DRAG, ENGINE, ENGWGT, PERFRM, TURN) 

1-3 

1 



ft 


ft 


ft 


ft 


ENGWGT (ENGINE, HOPWSZ, RCWSZ 
WGHT (LOAD, ENGSZ , ENGWGT, TAIL) 
OUTPUT (CLIFT) 


AEROUT 

(CLIFT, 

DRAG) 




PERFRM 

(ACCEL, 

CLIMB, 

DLAND, 

TAKOFF , 

.TAXI, TURN, X RANGE) 

RGBAL 

(AEROUT, 

CTAER, 

ENGSZ, 

ENGWGT, 

FLAPS, OUTPUT, PERFORM, 


SIZE, WGHT) 

GACOST (ASPEED, ENGINE) 

PNOYS (ASPEED, ENGINE, GEARBX, ZNENG) 

It may be noted, for example, that subroutine PERFRM is also called by 
ENGSZ and RGBAL, and that ENGSZ is called by WGHT. That is, there exists a 
very strong and complex connection between the various subroutines and the 
final effect of changing a parametric value is usually impossible to predict 
a priori. A detailed flow chart for the MAIN program is presented in Figure 
1.3.2. 
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FIGURE 1.3.1 - MAIN PROGRAM & SUBROUTINE STRUCTURE 
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I wgi=wc; 
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HPMSL1 =HPMSLS 
SWSL 1=SWSLS 
DPR0P1=DPR0P 
KIVRIT1=KKRITE 
CKVT1=CKVT 
CKHT1=CKHT 
CKF1=CKF 
CKBM1=CKBM 
CKN1=CKN 
CKTP1=CKTP 
CKW1=CKW 
CDNAC=0. 
GRFE1=GRFE 
DC LMT 1 = DC LMTE 
DELTF.1=D0LTE 
DCDOI 1=I)CD0TE 
WCFLP1-WCFLAP 
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DI’LPLI), DCLMLE , WCIT.AP , (TOC , 
DELLEO, BENGOB, BTEOB , DCLMTli , 
DCDOTE , PEI.TEO 


WRITE ((>)KWCI), 
ACLS ,AC1)CDR 
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1.3.2 Subroutine BIV - Linear Interpolation 
in Two Independent Variables 

A utility routine performing a linear interpolation in stored data 
of the form 


Z., = Z.. (X., Y. ) i = 1, 2, . . ., N 
3k 3k 3 k i 

j = 1, 2, . . ., 

Interpolation only is permitted. If an independent variable value falls 
outside the stored range, an error exit is made and the independent variable 
values being employed together with their boundary values are printed out. 
Figure 1.3.3 presents a detailed flow chart for this subroutine. 
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FIGURE 1.3.3 SUBROUTINE BIV 
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1.3.3 Subroutine INTS - Double Precision 


Finite Difference Integrator 
The calling sequence is 

CALL INTS (T, M, L, E, B, C, HMA, HMI, BET, DERIV) 

This utility routine is a finite difference integrator, performed in double 
precision, of a system of M simultaneous first-order differential equations 
which are defined in external subroutine DERIV. The non-zero components of 
T (100) are related to the state variables in DERIV. The other parameters 
in the calling sequence are input and are associated with the numerical aspects 
of integration (error magnitudes, step sizes, etc.). 

A detailed flow chart for INTS is provided in Figure 1.3.4. 
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FIGURE 1.3.4 SUBROUTINE INTS 








FIGURE 1.3.4 SUBROUTINE INTS 
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j ENTRY 

INTM \ 

1 

t 


c k -» y -r 


0 


X0=T(2) 

XX(2)=X0+H*0.5D0 

XX(1)=XX(2) 

! XX(4)=X0+H 
ixX(3)=XXf4) 


-►DO 16 JJ=1,4 


T(2) =XX (JJ) 1 

~T^ 

► DO 12 1=1, N 


* I 


11=3+1 
IJ=I I+N 
IK=IJ+N 
IL=IK+N 
AK=H*T(IJ) 


c 


JJ 


11 2 


T 


> 


T(.IK)=T(II) 

T (IL)=AK 

T(II)=T(IK)+0. 5DO*AK 


HI 


CONTINUE 1 


- jT(II)=T(IK)+T(IL)+AK)/6.0DO 


■ |t(II)=T(IK)+AK | - 


|T(II)=T(I 

K) +0. 5D0*AK | 


4 

f 

[T(IL)=T(IL)+2.0D0*AK j 


ICALL DliRIV 


ICR 




- { 16 CONTINUE I 


■0 




-►RETURN 
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DO 23 1=1, N 


11=3+1 

IJ=II+N5 

IK=IJ+N2 

IL=IK+N2 

IM=IL+N2 

IN=IL+N 

IP=II+N3 


11=3+1 

IJ=II+N2*(J+2) 

T(IJ)=T(II) 


21 CONTINUE 


RETURN 


T(IP)=19.0D0*T(IM)-5.0D0*T(IL)+T(IK) “ 
T(II)=T(IN)+H/24.D0*C55.D0*T(1M)-59.D0»T(IL)+37.D0*T(IK)-9.D0*T(IJ)1 

— I ~ ” 

23 CONTINUE I 


XS(2)=T(2) 
X0=XS (2) 

T (2)=X0+H 


CALL DERIV 


‘DO 24 1=1, N 


WRITE (6) 


11=3+1 
IJ=II+N2 
T(IJ)=T(II) 

IL=II+2*N5 

IK=II+N 

IP=II+N3 

T(II)=T(IL) +H/24 . DO* (9. ODO*T (IK) +T(IP) ) 


24 CONTINUE 











FIGURE 1.3.4 SUBROUTINE INTS 
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DO 27 1=1, N 











FIGURE 1.3.4 SUBROUTINE I NTS 


5 



RETURN 
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1-3.4 Subroutine ITRLN - Linear 


Interpolation in One Independent Variable 
This routine performs a linear interpolation in stored data of the 

form 

Y. - Y. (X. ) i'=l# 2, . . . , N 

1 li 

The calling sequence is 

CALL ITRLN (AX, AY, X, Y, N) 

This subroutine returns a value for Y corresponding to an input 
quantity X. The input parameters are the N pairs AX (I) and AY (I), 
and AXOIP must increase monotonically. If X is less than AX(1) or 
greater than AX(N), the subroutine extrapolates for Y(X). 

A deailed flow chart for ITRLN is presented in Figure 1.3.5. 
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1.3.5 Subroutine ITRMHW - Location of Root 
by Newton- Raphson Method 

This utility routine finds a zero of the function 

E = f(D) 

The calling sequence is 

CALL ITRMHW (ERROR, ERRM1, DRIVER, F, FF, JC, JX) 

This subroutine determines a zero to a function defined externally . 
Inputs are ERROR, the current (non-zero) value of the dependent variable; 
DRIVER, the current value of the independent variable; and F, a multiplier 
near unity. Outputs are ERRMl and DRIVER, the augmented values of the 
dependent and independent variables, and JC, the counter. FF and JX 
are not used. 

A detailed flow chart for ITRMHW is presented in Figure 1.3.6. 



FIGURE 1.3.6 SUBROUTINE ITRMHW 
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RETURN 






» 


i 


» 


I 


» 
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1.3.6 Subroutine MAXMHW - Maximum of a Function 
of One Independent Variable 

This utility routine determines a local maximum of the function 

Y * f(D) 


The calling sequence is 

CALL MAXMHW (PARAM, PRMMl, DRIVER, F, FF, KC, KX) 

This subroutine determines the maximum of an input function Y (X) "PARAM (DRIVER, 
which is defined externally. F and FF are input multipliers near unity, 
and KC is an output interaction counter, while KX changes from 0 to 1 when 
the maximum is determined. The previous value of Y(X) is PRMMl, and 
DRIVER is both input and output value of X. 

A de tai led flow chart for MAXMHW is presented in Figure 1.3.7. 
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FIGURE 1.3.7 SUBROUTINE MAXMHW 




1.3.7 Subroutine OUTPUT - Program Print Output Routine 


This routine provides print output of the aircraft characteristics. 
Calling sequence is 


CALL OUTPUT 

This subroutine begins with thirteen common block statements, and it 
includes 34 FORMAT statements. The subroutine is called by MAIN for the 
purpose of printing over 100 input and output figures related to geometry, 
weights, aerodynamics or the aircraft design. 

A detailed flow chart for subroutine OUTPUT is presented in Figure 1.3.8. 





DBARN=ANC/3. 14/XLN 

ELN=XLN 

SNX=ANAC*ENP 


VFTX=VFTP/2. 


WRITE (6) VDMIN, 
VMO,EMMO,ULF,EMLF, 
GLF , WEP , WPE I , WFSS , 
WPROP, WP,WW, WHT, 
WVT 


WRITE (6) 
WLG.WPES 


SfC 


tOUTINE OUTPUT 








FIGURE 1.3.8 SUBROUTINE OUTPUT 













RETURN 








1.3.8 Subroutine TPALT - Atmospheric Properties Routine 


This routine provides characteristics as a function of altitude. The 
calling sequence is 

CALL TPALT (ALTZ , ALT, PO, FKALT, TO GO, XKV) 

This subroutine relates static pressure, temperature and gravity, and kine- 
matic viscosity in ft 2 /sec, (PO, TO, GO, XKV) to the altitude. ALTZ is geo- 
metric altitude, ft., and ALT is potential altitude, ft., while PO is measured 
in lb per sq in., TO in deg R, and GO in ft per sec per sec. 

If PO is in input, ALTZ and ALT are output, and vice versa. FKALT determines 
whether geometric or geopotential altitude is used. 

A detailed flow chart for TPALT is presented in Figure 1.3.9. 



FIGURE 1.3.9 SUBROUTINE TP ALT 


Subroutine TPALT (ALTZ, ALT, PO, FKALT,T0,G0, XKV) 
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FIGURE 1.3.9 SUBROUTINE TP ALT 


2 








FIGURE 1.3.9 SUBROUTINE TP ALT 


3 


c 


^6p) 


T 




— ►DO 480 IX=1 , 8 

I 

( ALTB(IX-t-l) 

y 

^80 CONTINUE | 


TljA8K=TB(iXJ 
TSLOPE (IX) 
OALT=A-ALTB(IX) 
T=TBASE+TSLOPE*DALTl 


TSLOPE 


> 


|P=PB ( IX) *EXP (- . 0187434*DALT/TBASE)| 


|p=PB(IX) * (TBASE/T) ** (. 018743424/TSLQPE)! 


XKV= . 270558E- 06*T*SQRT (T) / (i . +198 . 72/T) ) 
G0=32. 1741+A7* (-309408E-05+AZ 
T0=T 
P0=P 
|ALTZ=AZ 
|ALT=A 


RETURN 
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1.3.9 Subroutine BILINE - 




Lin ear Interpolation, One Independent Variable 

This is a utility routine performing linear interpolation in stored 
data of the form 


* Z^(X^) s i * 1, 2, . • • N 


A det ai led flow chart for this routine is provided in Figure 1.3.10. 




» 


» 


1-3 

34 


» 




G'S' 


I 



BILINE 



( 20 (^ 


NY 


c 



f 

|JY=JX+NX | 

y 

t 


>T 


L =i 
JX=YI 
|J1=J2+1 
|j2=Jl+NY-l 


TN 


I 


* > z= °-° h 


jz=Cl* (T (JY+1) -T(JY) )+T (JY) | 

▼ 

RETURN 


f 


'DO 220 J*l,4 

1 


I Z=Z+C (J)*T (JY) 
JY=JY+1 


220 CONTINUE | ~ 


-►RETURN 


FIGURE 1.3.10 - SUBROUTINE BI LINE 
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biline 


2 



[K=K+3 

[kxJ 


r 


t ™ y j 

y 


JY=JX1+NY* (JY-IPX-2) 

JX=JY+NY 

Z=T (JY) +C1* (T(JX) -T(JY) ) 

Z= (X-T (JX1) )/ (T (JX1+1) -T (JX1))* (T(JY+1)+C1* 
(T (JX+1) -T (JY+L) ) -Z)+Z 

▼ 

RETURN 


♦ * ~ 

tlY=J2+l+(JX-IDX-3)*NY+JXl-Jl | 

r 

DO 550 M=l,4 



| 550 CONTINUE I 




iz=0.0 I 

DO 700 J*l,4 
jz=Z+C(J)*Y(J)l 
RETURN 



[105] 


t 






BILINE 


3 

















BILINE 


C4=RB/P5*D2/P3*D3 

C1=RA/P1*D2/P4*D3 

C2=-RA/P1*D1/P2*D3+RB/P2*D3/P5*D4 

C3=RA/P2*D1/P4*D2-RB/P2*D2/P3*D4 






1.3.10 Subroutine BIQUAD - 
Quadratic Interpolation, One Independent Variable 


This is a utility routine performing quadratic interpolation in one 
independent variable using data stored in the form 

Z i “ z i (x i )j i - 1/ 2, . . . n 

A detailed flow chart for subroutine BIQUAD is presented in Figure 1.3.11. 
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BIQUAD 


2 



DO 1010 J=J1, J2 



DO 2020 J=l,4 


D(J)=X-XC (J) 


2020 CONTINUE 


C(l)= (RA/P(1))*(D(2)/P(4))*D(3) 

C(2)=(-RA/P(1))*(D(1)/PC2))*D(3)+(RB/P(2))*(D(3)/P(S))*D(4) 
C(3)= (RA/P(2))* (D(1)/P(4)*D(2)-(RB/P(2))*(D(2)/P(3))*D(4) 
C(4)=(RB/P(S))*(D2)/P(3))*D(3) 
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1.3.11 Subroutine MAXBND - 


Maximum Value of a Variable 

This is a utility routine which determines the maximum value of a 
variable in the interval DM1N to DMAX. A detailed flow chart for MAXBND 
is presented in Figure 1.3.12. 



MAXBND 
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MAXBND 


2 





KM=0 

KP=0 

KA=0 


\ 

RETURN 




I 


I 


> 
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1.3.12 Subroutine UNINT 


Four Point Smooth Interpolation 

This is a utility routine which performs a smooth four point interpo- 
lation in stored tabular data of the form 

U i " Y i (X i )l 1 ■ 2, . . . N 

A detailed flow chart for UNINT is provided in Figure 1.3.13. 
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UNINT 




RETURN 









